HEAT ASYMPTOTICS WITH SPECTRAL BOUNDARY 

CONDITIONS 
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Abstract. Let P be an operator of Dirac type on a compact Riemannian 
manifold with smooth boundary. We impose spectral boundary conditions 
and study the asymptotics of the heat trace of the associated operator of 
Laplace type. 



1. Introduction 

Let M be a compact Riemannian manifold of dimension m > 3 with smooth 
boundary dM. Let Ei be unitary bundles over M and let 

(1) P:C°^{Ei)^C°^{E2) 

where P is first order partial differential operator. Let P* be the formal adjoint 
of P. We say display is an elliptic complex of Dirac type if the associated 
second order operators 

Di := P*P and D2 := PP* 

on C^{Ei) and on C°°{E2) are of Laplace type - i.e. if these operators have 
scalar leading symbol given by the metric tensor. If Ei = E2 and if P = P*, 
then P is said to be an operator of Dirac type; it is convenient, however, to 
work in this slightly more general context. 

We impose spectral boundary conditions; these were first introduced by 
Atiyah, Patodi, and Singer in their study of the index theorem for manifolds 
with boundary. Let 7 be the leading symbol of the operator P. Then 7 + 7* 
defines a unitary Clifford module structure on Ei ® E2. We choose a unitary 
connection V = Vi © V2 on Pi © E2 so that 

(2) V(7 + 7*) = and (Vs, s) + (s, Vs) = rf(s, 5); 
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such a connection always exists and is said to be a compatible unitary 
connection. We note that equation (H) does not determine V uniquely; there 
are many compatible unitary connections. 

Use the inward geodesic flow to identify a neighborhood of the boundary with 
the collar dM x [0, e) for some e > 0; if ?/ G dM, then the curves y(t) = {y, t) 
are unit speed geodesies perpendicular to the boundary. If ?/ = {y^, ...,y"^~^) 
are local coordinates on dM, let x = {y, x"^) be local coordinates on the collar 
where x"^ is the geodesic distance to the boundary. Let := dm is 
the inward geodesic normal vector field on the collar. Let be covariant 
differentiation with respect to d^. We may decompose 

l<fj,<Tn 

where ip is a 0*^ order operator. Since we do not assume that the structures are 
product near the boundary, the connection 1 form of V, the leading symbol 
7, and the endomorphism ip can depend on the normal variable. Relative to 
a local frame on the collar which is parallel along the normal geodesic rays 
which are perpendicular to the boundary, we have Vm = dm- We set x"^ = 
to define a tangential operator 

Let be an auxiliary self-adjoint endomorphism of E\\qm- Let 

A:=\ [b, + 5*^ + on (^^4,,) . 

Here, the adjoint of -Bq is taken with respect to the structures on the boundary. 
The operator A is a self-adjoint operator of Dirac type on C°^{E^qm)- Let the 
boundary operator i3, which we will use to define the boundary conditions for 
the operator P, be orthogonal projection on the span of the eigenspaces for the 
non-negative spectrum of A. Denote the realization of P and the associated 
self-adjoint operator of Laplace type by 

Prs and Di,^ := {PbYPb- 

Let F e C°°(i5i) be an auxiliary function we use to localize the problem; F 
is called the smearing function. Results of Grubb and Seeley [0, [T^, |TB| show 
that there is an asymptotic series as t J, of the form: 

(3) Tri2(Fe-*^i'«) ~ ^ ak{F, Di,B)t^''-"'^^^ + 0{t-^/''). 

0<fc<m-l 

(There is a complete asymptotic series with log terms, but we shall not need 
this fact as we shall only be interested in the first few terms in the series). 
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The coefficients in equation are locally computable. We shall deter- 
mine ao, ai, and 02 in this paper. Our purpose is at least partly expository, 
so we will present several different techniques which yield information about 
these asymptotic coefficients. 

We adopt the following notational conventions. Roman indices i and j will 
range from 1 to m and index a local orthonormal frame for TM; Greek indices 
will index a local coordinate frame. Near the boundary, we choose the frame 
so that Cm is the inward unit geodesic normal vector; we let indices a and b 
range from 1 through m — 1 and index the corresponding frame for the tangent 
bundle of the boundary. We adopt the Einstein convention and sum over 
repeated indices. We let ';' denote multiple covariant differentiation of the 
tensors involved. Decompose 

D, = -{g''''d,d, + a^d, + b). 

Let r be the Christoffel symbols of the Levi-Civita connection on M. There 
is a canonical connection on the bundle Ei and there is a canonical endo- 
morphism E of the bundle Ei so that 

Di = -{Tr(^V') + F}; 

we refer to [jl4| for details. If Ei = E2 and if P = P*, then is a unitary 



connection; however need not vanish in general so is not in general 

a compatible connection. Let u be the connection 1 form of ^V. We have 



1 

(4) E:=b-g''^{d,u^ + uj,uj^-u^T,^''). 



c^5:=7T^?-5K + ^7'^"r^/) and 



Decompose P = 7iVj + ip and let ip := 'JrJ'ip- Let 71 be the scalar curvature 
of the metric on M. The main result of this paper is the following: 



Theorem 1.1. We have 

1. ao(P,Z^i,S) = (4vr)-™/2/MTr(P). 



2. ai(P,Di,i3) 



r(f) 
r(i)r(i2±i) 



1 



{An 



,-(m-l)/2 



/aMTr(P). 



2{m-2) 



r(|)r(f ) 
2r(i2i±i) 



F + 



3r(|)r(f ) " 



LaaF 



;m 
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U k < m, then there exist locally computable invariants a^^ and af*^ so that 
a,iF,D,,B) = (47r)-™/2 /" Tr (af (F, Di, x)) 



JdM 



We have included a normalizing factor of (47?)^™/^ to simplify the formulas 
for the local invariants and af*^. We use dimensional analysis to see that 
the invariants a^^ are homogeneous of weight k and the invariants af^ are 
homogeneous of weight — 1 in the jets of the symbols of P and P*. The 
invariants a2j+i vanish. We use Theorem 4.1.6 0] to see that: 



(5) (F, Di,x) = Tr {F} and af (F, a;) = Tr {F(-7^ + F)}. 

6 

The bundles Fi and F2 are distinct; we must use 7^ to identify Fi and F2 
near the boundary. This observation reduces the number of invariants which 
are homogeneous of weight 1; for example, Tr (tp) is not invariantly defined. 
There are universal constants so that 

af\F,Di,B) = h{m)Ti{F) 

af ^(F, Fi, B, y) = Co{m)F{^ + r) + c,{m)F{^ - r) 

+C2{m)FQ + Cz{m)FLaa + Ci{m)F.rn- 

In contrast to the situation when the boundary operator B is local, the con- 
stants exhibit non-trivial dependence upon the dimension. 

We will use three different methodologies to compute the unknown coeffi- 
cients. In Section |^, we use results of Grubb and Seeley [|1^ for structures that 
are product near the boundary to compute the constants 61 (m) and co(m); see 
Lemma for details. In Section S, we use functorial properties of these in- 



variants to determine the coefficients co(m), ci(m), and C2(m); see Lemma 
for details. In Section ^, we use computations on the ball to determine the 
coefficients &i(m), c^{m), and c^{m)\ see Lemma ^?l|for details. As a check on 



our methods, we give two different derivations of the relation 



(6) (m — 2)c4(m) + (m — l)c^{m) - 

6 

in Sections ^ and ^. Our purpose in this paper is partly pedagogical; we wish to 
illustrate different methodologies which can be used to study these invariants. 



2. Product Formulas 
We use results of Grubb and Seeley [|18[ to show: 
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Lemma 2.1. 

1. We have bi{m) 

2. We have Co{m) 



r(f) 



r(i)r(i2±i) 



(47r)-(--i)/2. 



Proof of Lemma |2.1| : Suppose that P = 7m (Vm + B) where B is a self- 
adjoint tangential operator of Dirac type with coefficients which are indepen- 
dent of the normal variable; we take Q = so A = B. In this setting, we 
say the structures are product near the boundary. Let Di be the associated 
operator of Laplace type on the double. We ignore the effect of the spectrum 
and define: 



r]is,B) :=Tr^.(i?(i?2)— 1) 
as,D,) ■.= Ttl^{D~^) 



C(2s,i?) :=Tri.((i?2)-) 



We refer to Theorem 2.1 [|l^ for the proof that: 



4^r(i)r(s + i) 



(7) 



Tis+W 



4r(i)r(s + i 



-ri{2s,B) 



where the remainder R is regular away from s = 0. Expand 



k 



The invariants a^^ vanish. Formulas of Branson and Gilkey show that 
(8) a;'(5) = -(47r)"(""^)/2^m-2) /" Tr(V'). 

JdM 

One can use the Mellin transformation to relate the asymptotics of the heat 
equation to the pole structure of the eta and zeta functions. Let be a 
manifold of dimension n, let D^v be an operator of Laplace type on A^, and let 
Qat be an operator of Dirac type on A^. If the boundary of is non-empty, 
impose spectral boundary conditions. We then have, see for example Theorem 
1.12.2 0, 



(9) 



ak{DN) 



Res„_n-kT(s)C(s, Dn) and 



Res„ 



2 

_n — k 



..T{s)ri{2s-l,QN). 
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The following identities now follow from equations (|^) and @: 

(10) a„(l, D,, B) = ^a„(l, A) - ^— t-<_i(1, 5) if n = mod 2, 

2 2[m — n)T[^) 

(11) a„(l, Z^i, i3) = - _ l)a„_i(l, 5^) if ^ ^ 1 niod 2. 

We prove the first assertion of Lemma ^]1] by using equations (||) and (|11]) to 
compute: 

(47r)— /2af'(l,A,S,l/) 



We prove the second assertion of Lemma ^Tl] by using equations (P) and (|10D 
to compute: 



(47r)~™/2af^(l,Z}i,S,y) 

= ^ ^ , r47r)-('"-^)/^a?(l,E,l/) 

2(m-2)r(i)^ ^ n > .y; 

= (47r)-™/2Tr (^/>). □ 



3. FuNCTORiAL Method 



The invariants have many functorial properties. We use these properties 
to establish the following result. 



Lemma 3.1. 



1. We have C2{m) = 0. 

2. We have co(m) = |. 

3. We have ci(m) = 0. 

4. We have (m — 2)cJm) + (m — 1)03(171) = — 



Proof of Lemma |3.1| : We take Ei = E2 and let P be a formally self-adjoint 
operator of Dirac type on C°°{Ei). Let Q{e) := + e where e is a real 
parameter. For generic values of e, keT{A{e)) is trivial. For such a value, the 
boundary conditions determined by the boundary operator B{e) are locally 
constant and thus 02 is independent of e. The first assertion of Lemma 
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now follows; this implies that a2(-) is not sensitive to the particular boundary 
condition chosen among the family we are considering. 

For the remainder of the proof of Lemma let T := S"^ x ... x S*^ be the 
m — 1 dimensional torus and let M := T x [0, 1]; we give T and M the canonical 
product flat metrics and let {xi, ...,Xm) be the usual parameters. We identify 
T with T X {0} in M. Let L be the line bundle defined by a non-trivial Z2 
valued representation of the fundamental group of T and let V be the trivial 
bundle of dimension 2™ with coefficients in L. Let 7^ be real skew-adjoint 
matrices acting on V which satisfy the Clifford commutation relations 

lilj + Ijli = -26ij. 

The twisting defined by L ensures that the kernel of the associated tangential 
operator B is trivial. Let Vj := di define a compatible unitary connection on 
C°°(y). Let Ao := i'jidi)'^ = —df be the associated operator of Laplace type 
on C°°{y). Let / = f{xm) be a smooth real valued function on M which 
vanishes identically near Xm = 1- 

We use the index theorem to evaluate the coefficient Co(m). Let 

Q := -fidi + f-fm, Q* ■■= lidi - fjm, 

:= Q*Q = A- 2-1„aafda + f- fym, 

D2 := gg* = A + 27^7,/9, + P + f,m. 
We use equation @ to compute: 

t^l.a = 7m7a/, ^l,m = 0, Ei = f-rn + (m - 2) f , 
^2,a = -Imlaf, ^2,ni = 0, E2 = - f-m + (m - 2) f . 

We have = f and ^2 = —f- The local formula for the index shows that the 
super trace vanishes for n ^ m. Since m > 3, we have 

= a2il,DuB)-a2il,D2,B) 

= dim(y){2 / f,^ + Aco{m) [ /} 
Jm Jt 

= dim(V)(-2 + 4co(m)) / /. 

Jt 

This shows that Co(m) = | and proves the second assertion of Lemma |3.1| ; 
note that this value agrees with the result obtained previously in the proof of 
Lemma |2.1j . 

Next we study the coefficient ci. Let 



Pf ■■= lidi + V^/77 
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this operator is formally self-adjoint. Note that Af = —•ymla^a so the bound- 
ary operator B is independent of the function /. We expand 

j2 



Since the 7^ were real, complex conjugation preserves the boundary conditions 
and intertwines Pj with P'^j- Thus Pj and P^j- are isospectral so 

(12) a,{l,PlB) = a2il,P^f,B). 

We use equation (^) to see the interior integrand vanishes by computing: 

oomif) = V^f and E{f) = 0. 

Since i'f = \/^f, we may use equation (|T^) to show ci = by computing: 

= a2(l,P;,i3)-a2(l,P!pi3) =4v^ci(m) / /. 



T 



We can also show ci(m) = using gauge invariance. Let P = e ^^^-^Poe^^-^ 
be defined by a global unitary change of gauge. We have 



Thus Tr {ijjf — ijjj) = — 2v— T/;m dim(V). As 02 is gauge invariant, ci(m) = 0. 

To prove the final assertion of Lemma |3.1| , we shall need the following tech- 
nical result; we postpone the proof until the end of this section. 



Sublemma 3.2. Let ds'^{e) = e^"^^ ds^ and let P{e) := e ^-^P. There exists a 
compatible family of unitary connections V(£:) so that 

i^{e) = e-^4m-'^e\f,a,. 



We use Sublemma ^]2| to complete the proof of Lemma |3TT| . Let Pq := 7iC?j 
and let 

ds\e) = e^'^ds'', Pie) := g-^M^/Poe-'^^'^^^^ 
dvol{e) = e'^'^dvol, P{ey = e^-P("^)~"^)ef p*(A^-<^{^))^f _ 

We fix the inner product on V. Let a{m) = and let /5(m) := 
Then a{m) + f3{m) = 1 and a{m) = f3{m) + m so the metric determined by 
the leading symbol of P{e) is ds'^{e) and P{e) is formally self-adjoint. We 
assume that / = f{xm) vanishes on the boundary of M and that / vanishes 
identically near Xm = 1- Thus the leading symbol of A{e) is independent of e 
so by defining 6(e) appropriately, we can assume the boundary operator B{e) 
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is independent of e. Let D{e) := -P(£)^. Let 6 := ^|£=o We compute the 
variation: 

5Tri2(e-*^(^)) = -tTi L'4{6D{e)}e~'^'') 

= -2m L'2{{SPie)}Poe-'^^^) = 2^X1^2 (/Aoe"*^") 

= -2t9iTr^2(/e-*^")- 
We equate coefficients in the asymptotic expansions to see 

(13) 6a,{l, D{e), B) = {m - 2)a,{f, Aq, B). 

If £ = 0, then D = Ao and af = 0. We are interested in the coefficient of 
£f;mm- Siucc u{0) = 0, siuce r(0) = 0, and since E{0) = 0, we may compute 

biD{e)) = -f3{m)ef,mm + 0{e^) 

a"'{D{e)) = (m - 2a{m) - 2f3{m))ef.m 

/^/ XX m — 2a(m) — 23(m) ^ , ■y^r ^ ^ 

uJm{D{e)) = ^-^ + cu^le'^^^ Aq) , 

E{D{e)) = b{e) - dmUJrnie) + i?(e"2-^Ao) + 0(8^) 
m — 2a(m) ^ ^, ^ , ^, 



2 

We use results of [0 to see that 

5E{e-^'^Ao) = -2fE + ^(m - 2)f,u, 
5n = -2/7^ - 2(m - l)f.ii, and 

This permits us to compute the variation of the interior integrah 
(14) ={-«(m) + ^}|^/^.. 

We have P{e) = e-=^(Po - P{m)e-fmf;m)- Since ^/'(Po) = 0, Lemma ^ (2) 
imphes that ifj^e''^^ Pq) = — ^^^^ee~^-^/;j7i. As / vanishes on dM, 

Si:{P{e)) = -I P{m) + ^— U,^ and 



mm 



(15) 5^af^(D(e),i3) = | 



-(77i - l)c3(77i) - /5(777) - } I f-m- 



Recall that a{m) + l3{m) = 1 and that 02 (Aq) = 0. We use equations (|T3|), 
(|T^) , and (|T5]) . The final assertion of Lemma |3.1| follows from the following 
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equation: 

TtX ~1~ 2 TtX — 1 

{m — 2)c4(m) = —a{m) H (m — l)c^{'m) — f3{m) . □ 

3 2 

Proof of Sublemma |3.2| : We follow the argument given in |jl2|. Let M be 

an arbitrary Riemannian manifold and let 7 be a Clifford module structure on 
a vector bundle over M. Let 9^ and dx^ be local coordinate frames for the 
tangent and cotangent bundles. We define: 

Let V(£) := V + 0. Since the original connection is compatible, we have 

= v^(7)(rfx^) = + + r^.V(o). 

Let S^{e) := V ^{e)^{e) and let 8"^ := Sfj^ldx"). We then have 

Fix xq G M and choose the local coordinates near xq so gfj,u{xo) = 5^^. Then 
Tf,y''{xQ) = ^g'^^id^g^.r + d^g^r - drg^u){xo) 



We must show £^{xq) = 0. If yU 7^ z/, we compute: 
t 1 

= e'""^ {-ef.^^y + ef-al,,5ya}{,xo) = 0. 
If /i = z/, then (don't sum over /i): 

^;(^o) = e-4-5/;^7^ + [e„r] + Y.(^,.'{e) - r,.^(0))7.)|(xo) 

This shows that V{e) is compatible. Let c = 2. Then 6'^ + 6^* = and V(e) is 
unitary. We complete the proof of Sublemma |3.2| by computing 



^(^) = e-^^l^ - ^£(2/^,7^7^7^ + cf,,Y) 
= e-'^!^tP - ^ef,^^>'{2m - 4 + c) | . □ 



heat asymptotics with spectral boundary conditions 
4. Computations on the Disk 
We perform computations on the disk to prove: 
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Lemma 4.1. 

1. We have 6i(m) 

2. We have c^^m) 

3. We have cJm) 



r(|)r(22±i) 



(4 



TT 



-(m-l)/2 



1 - 



3r(|)r(f ) 
4r(2i±i) 



m-i 1 _ r(^) 

2(m-2) [ 2 V2'' r(2i±i) 



4. We have (m — 2)c4(m) + (m — l)c3(m) 



6 



Proof: Let M be the unit ball in with the usual metric. If r G [0, 1] is 
the radial normal coordinate and if (iS^ is the usual metric on the unit sphere 
S"^~^, then ds"^ = dr"^ + r'^dT? . The inward unit normal on the boundary is 
—dr- The only nonvanishing components of the Christoffel symbols are 

Tafec = ~^abc and Tabm = ""^afe; 

the second fundamental form is given by Tabm = Lab- Here Tabc are the 
Christoffel symbols associated with the metric dT? on the sphere 5"™"^ and 
tilde will always refer to this metric. 

The spin representation 7 is an irreducible representation of the Clifford 
algebra; we refer to for details. Let P = "^^dy be the Dirac operator on the 
ball; we take the flat connection V and set ip = We suppose m even (there 
is a corresponding decomposition for m odd) to find a local decomposition: 






"1 ■ lm-1 



-1 ■ 7(Vi) 




and 



-1 ■ 1 



m— 1 



We stress that 7^^-) are the 7-matrices projected along some vielbein system. 

Decompose Vj = Cj + Uj where ojj = jTjki'jk'Ji is the connection 1 form of the 
spin connection. Note that 











+ 7^Sabi-l}m)lim), 
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Let P the Dirac operator on the sphere. We have: 

d m-l\ ^ \ f ^^^P 



dx^ 2r J '^'^> r \ -V^P 

Let ds be the dimension of the spin bundle on the disk; dg = 2™/^ if m is even. 
The spinor modes on the sphere are discussed in We have 

PZ^^\n) = ±(^n + ^^^^ Z^^\n) for n = 0, 1, 

^ / \ A- -7H ( m + n- 2 

dnim) := dimZj. = -d ' 



2 V n 



Let Ju{z) be the Bessel functions. These satisfy the differential equation 
and functional relations 111 511: 



dz^ z dz \ z^ J 

z-^JJz) + pJJz) = zJy-liz), and 
dz 

z—J^{z) - vJy{z) = -zJ^+i{z). 
dz 

Let P(p± = ±A(y9-|- be an eigen function of P. Modulo a suitable radial nor- 
malizing constant C, we may express: 

Let ^ll^) := -l'(:n)lU and let ^V, := V„ - \Lab^l\^y Then is a com- 
patible unitary connection for the induced Clifford modules structure ^7; see 
p!3| for details. We may express the tangential operator B in the form: 

^ = -7em)7(m)Va= 7em)Va= 7(m)( ^a + -Labl(m)\ 
_P_n^ 

^- ^ 

Thus in particular B = B* . We take 6 = 1^- We then have: 



A 



2 

-P 
P 
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The eigenstates and eigenvalues of A then are given by: 

^ I ^t„) ) = - (n + "^^^^—^ 1 I 1 and 

^(n) 




^ I ^[n) 1 = 1 ^ + 1 I ^7n) 1 forn = 0, 1, 



The boundary condition suppresses the non-negative spectrum of A. We use 
equation (|I7p to see that the non-negative modes of A are associated with the 
radial factor J„_|_^_i(Ar). Hence the implicit eigenvalue equation is 

TfX 

(18) Jp(A) = where p = n+ — -l. 

The first and third authors developed a method for calculating the associated 
heat-kernel coefficients for smearing function F = 1 in ^, |10]; they gener- 
alized this method to deal with F = F{r) in ||lT|. We summarize the essential 



results from these papers briefly; in principal one could calculate any number 
of coefficients. We first suppose that F = 1. Instead of looking directly at the 
heat-kernel we will consider the zeta-function ([s) of the operator and use 
the relationship provided by equation (^ between the pole structure of the 
zeta function and the asymptotics of the heat equation: 

(19) Ofc = Res „^ m-k T(s)C(s). 



2 



Thus to compute ao, ai, and 02, we must determine the residues of the zeta- 
function ({s) at the values s = s = and s = ^ — 1. We use the 
eigenvalue equation (ITB|) to express 



(20) C(^) = 4 d^{m) / ^A;-^^- In J,(fc), 

where the contour C runs counterclockwise and encloses all the solutions of 



T8|) which lie on the positive real axis. The factor of four comes from the four 
types of solutions in (p!6D and (|^). As it stands, equation (pO|) is well defined 
only for JRs > m/2, so the first task is to construct the analytical continuation 
to the left. We define a modified zeta function 



'C 

the additional factor has been introduced to avoid contributions coming 
from the origin. Since no additional pole is enclosed, the integral is unchanged. 

The behaviour of C''"^^) as n — 00 controls the convergence of the sum over 
n; different orders in n can be studied by shifting the contour to the imaginary 
axis and by using the uniform asymptotic expansion of the resulting Bessel 
function Ip{k). To ensure that the resulting expression converges for some 
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range of s when shifting the contour to the imaginary axis, we add a small 
positive constant to the eigenvalues. For s in the strip 1/2 < 3fJs < 1, we have: 



c 



sin(7rs) 



TT 



dk{e - e^'— In k-PIp{k). 



We introduce some additional notation dealing with the uniform asymptotic 
expansion of the Bessel function. For p ^ oo with z = k/p fixed, we use results 
of iTH to see that: 



f211 



lp{zp) 



1 



(1 + ^2)1/4 



pi 

1=1 ^ 



where 



t = 1/Vl + z^ and r] = Vl + + \n[z/{l + Vl + z% 

Let uo{t) = 1. We use the recursion relationship given in to determine the 
polynomials Ui{t) which appear in equation (p]) : 

WmW = lt\l - e)u[{t) + 1 ^ rfr(l - 5t')ui{t). 
We also need the coefficients Dmif) defined by the cumulant expansion: 



(22) 



In 



E 

1=1 



ui{t) 



pi 



E 

9=1 



pq 



The eigenvalue multiplicities dn{m) are 0{n"^~'^) as n ^ cx). Consequently, 
the leading behaviour of every term is on the order of |j-2s-g+m— 2. ^^j^^g on the 
half plane 3?s > (m — 3)/2, only the value q = I contributes to the residues 
of the zeta-function. We have Di{t) = |i — -^t^- We use equation ( ^1]) to 
decompose 



A'llis) + A^^\s) + + R^''\s), where 



A^r\s) 



sm ITS 

TT 

sinvrs 

TT 

sinvrs 

TT 



d 



e/p 



dz[{zpy -e^]-'—\n {z-PeP'^) 



dz[{zp) 



2 g2l-s 



dz 
d 



e/p 



In (1 + ^2) 



-1/4 



dz[{zpy 



dz 



e/p 



dz 



p 



The remainder R^^\s) is such that 'Yl'^=o^n{fn)R^^'^{s) is analytic on the half 
plane 3?s > (m - 3)/2. 
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Let 2F1 be the hypergeometric function. We have 

r(c) 



15 



2Fi{a,b; c; z) 



T{h)V{c-h) 



[ dtt^~\l - tY~^-\l - tz)-\ and 



dz [{zp) 



lV{s + {)V{l- s) 



e/p 



^ dz' 2 r(i + i) +^ 



l\ 2 I 2]-s-l/2 



We use the first identity to study A^^l{s) and ylQ"''(s); we use the second identity 
to study A^^\s). This shows that 



,-2s+ir(.-i) 1 _i 1 
2r(i) r(s) 2'^ 2' 2' 

1 1 

5 1 
"24f(s) 



e 



2 



We take the hmit as e — 0; the resulting zeta-function which appears is 
connected to the spectrum on the sphere. We define the base zeta-function 
and the Barnes zeta-function |^ 

00 00 
Csd{s) = 4:'^dn{m)p~^^ and (3(3, a) = ^dn{m){n + a)'" . 

n=0 n=0 

We then have the relation (s^'i^) = 2(isCB (2s, y — l). For i = —1, i = 0, and 
i = I, let Ai{s) = 4^^Q(i„(m)y4."''(s). We take the limit as e — > to see that 



(23) 



A_,is) 



4r(|)r(s + i) 



(24) Ao{s) = --Cs4s), and 



(25) 



AAs) 



Vis 



1 ^ / ^ 



2' [8r( 



i2r( 



We use the Mellin-Barnes integral representation of the hypergeometric func- 

V{a + t)V{b + t)V{-t) 



tions [0 to calculate A_i{s): 

r(c) 1 



2Fi(a,6; c; 2;) 



r(a)r(6) 27ri 



-(-.)*. 



r(c + t) 

We choose the contour of integration so that the poles oiV{a+t)V{h+t) /r(c+t) 
lie to the left of the contour and so that the poles of r(— t) lie to the right of 
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the contour. We stress that before interchanging the sum and the integral, we 
must shift the contour C over the pole at t = 1/2 to the left; this cancels the 
term — le"^'' appearing in the expression for A^i appearing above. 

This reduces the analysis of the zeta function on the ball to analysis of a zeta 
function on the boundary. We compute the residues of ({s) from the residues 
of Cb{s, a). Let d := m — 1. To compute these residues, we first express Cb{s, a) 
contour integral. Let C be the Hankel contour. 

CB(s,a) = ( M (^ + «)"' = (a + mi + ...+mrf)-" 



rn-s) r , e 



-at 



271 Jc (l-e-^y 

The residues of Cb{s, cl) are intimately connected with the generalized Bernoulli 
polynomials |I9| , 

-at °° / ±\n-d 

(26) 7r^ = (-i)'5:^^'!.'"w- 

^ ' n=0 

We use the residue theorem to see that 

(27) Res _Cb(^, a) = 

for z = 1, ...,d. The needed leading poles are 

1 



Res s=dCB{s, a) 
Res s=d-iCB{s, a) 
Res s=d-2CB{s,a) 



{d-l)V 
d-2a 



2{d-2)V 

\2o? -d - \2ad + ?>d^ 



24(c/-3)! 

We may now determine the residues of At s = only v4_i(s) contributes. 



We use equation (|23|) to see 

r(^) 1 



Res s=sC(s) = d. 



'4y?r(f + l)r(m-l)" 

We use the 'doubling formula' ~ ^ function, we use 

equation (0), and we use the observation Tr (1) = ds to check that we have 
derived the correct formula for oq: 



ao 



^ = (47r)-"^/2 / Tr (1) 

2-r (f + 1) ^ ' Jb^ 
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Next, we study the pole at s = This time, A_i{s) and Ao{s) contribute. 
We compute 



Res 
Res 



m — 1 
" 2 



2 



A_As) 



An(s) 



1 



2mr(!Ii±l)r(H^)' 



This imphes that 
1 



2™r(f ^ 



4r(m - !)■ 



r(f) 



4 Lr(|)r(H^) 



(47r) 



Tr(l). 



This determines bi{m) and proves the first assertion of Lemma [4.1| . Next we 
compute a2 and thereby determine Cs{m): 



Res „_m-2 74_i (g) 
Res 5=iri^Ao('5) = 

Res ,=inz:2Ai(s) = 



1 4 — m 
6 2" 



1 



m-2 \ ' 



8r(m-2; 

1 8 — 5m 

Q 2m 



m-2^ 



We can now prove the second assertion of Lemma [4.1| by determining C3(m) 
We compute: 



2 (m 



a2 



1 (^-i)r(i) 



3 2'"r(f ^ 



2 2"'T( ''"^"^ ' 



ill 



3r(|)r(f; 
4r(^) 



(4vr: 



|-m/2 



Tr(L„ 



The pattern of the universal constants involving a dimensionless constant and 
a combination of F-functions is evident. This pattern holds for the higher 



coefficients HIO 



We introduce the weighting (or smearing) function F = 1 — to determine 
C4(m); since E and TZ vanish and since F vanishes on the boundary, only the 
term involving F-m = 2 survives. (Note: we checked these computations by 
studying a more general function of the form F{r) = /o + fir^ + /2'"^ but 
omit details in the interests of brevity). We note that the radial normalization 
constant is given by C = l/Jp+i(A). We denote the normalized Bessel function 
by 

Jfc(Ar) := Jfc(Ar)/Jp+i(A). 
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Instead of the zeta function we consider now the smeared analogue: 

1 



(28) 



x)(p*{x)(p{x) 



X2s- 



Since F depends only on the normal variable, the integral in equation 
over the sphere S"^~^ behaves as in the case F = 1 so that 



oo „ 

4Vd„(m) / 



n=0 
1 



dk 
27ri' 



k 



-2s 



d 



drF{r)r{J^_^,{kr) + J^{kr))—\n J,{k). 



We set F{r) = 1 — and compute the radial integrals: 



K(Ar) + J,V(A^ 



2p2 + 3p + 1 1 



3A2 



so 



C(r2;.)=4f^rf.(m) / ^k^ 



2s 



2p2 + 3p + i 1' 
3P +3 



^lnJ,(A:). 



We use equation (^) to evaluate this expression; the second term is given by 
equation (^) and simple substitutions suffice to evaluate the remaining parts. 
The factor is absorbed by using s + 1 instead of s in equations p3|), (plf), 
and (pSf). The powers of p lower the argument of the base zeta-function by 1, 
by I and by 0. It is now a straightforward matter to compute: 



1 ^ 

4r(i)r(s + i)''^'^' 2' 



1 2s - i 

- H 

3 3s + 1 



4r(i)r(s + 2) 



■^Cs4s) - ^Cs4s + -) + .. 



2 1 

"3r(s + i)^^'^'+2^ 



sr 



1 ^/ 3, 



5 5, 



5 3, 

T2ra)^^'^2^ 



+ ... 



The coefficient oq gives the leading term in the expansion of the heat trace; 
this is correctly reproduced by the first term in A_i{r'^;s). We also confirm 
the invariant form of the next coefficient ai. We determine the value of C4(m) 
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by considering the pole at the value s = Note that F.„^ = 2. We compute 

1 m 



Res s=iB::i2r(s)y4_i(l — r ; s) 



12m -2 



(An) 



Res „_m-2 
*~ 2 



4 '2'm-2r((m + l)/2) 



(47r) 



-m/2 



Res „_ m-2 
*~ 2 



T{s)Ai{l -r^;s 



1 5m — 6 



(47r 



-m/2 



Tr(2) 



-m/2 



Tr(2) 



Tr(2). 



12 m - 2 

We sum the contributions to prove the third assertion of Lemma |4.1| by eval- 
uating C4. The final assertion of Lemma is now immediate. □ 



5. Conclusion 

In Section ||, we used results of Grubb and Seeley to study the setting when 
the structures are product near the boundary. In this setting, ip was self- 
adjoint and only the terms Tr (1) and Tr {tp) were non-trivial. We used these 
results to determine the coefficients 61 (m) and co(m). 

In Section |, we used functorial methods to determine the coefficients co(m), 
ci(m), and C2(m); the value we computed for co(m) agreed with that deter- 
mined in Section ^ We used variational methods to obtain a non-trivial linear 
relationship between the coefficients Cs{m) and c^^m). 

In Section |, we determined 61 (m), c^^m) and c^^m) by studying the ball in 
Euclidean space. The value of 61 (m) computed in this fashion agreed with the 
value which was determined in Section |^. The values for the coefficients C3(m) 
and C4(m) satisfied the linear relationship derived in Section ^. 

In Section we chose c = 2 to ensure that the 1 parameter family V(e) of 
connections was unitary; the second author originally missed this point and 
his error in choosing a different value of c led to a seeming contradiction, which 
has now been cleared up, between the methods of Section |^ and the methods 
of Section ^. We find it interesting that each of the three methods we have 
discussed gives some, but not all, of the coefficients and that none of the results 
of the three sections is a proper subset of the other. 
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